Calculation of the spectrum of the superheavy element Z = 120 
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High-precision calculations of the energy levels of the superheavy element Z = 120 are presented. 
The relativistic Hartree-Fock and configuration interaction techniques are employed. The correla- 
tions between core and valence electrons are treated by means of the correlation potential method 
and many-body perturbation theory. Similar calculations for barium and radium are used to gauge 
the accuracy of the calculations and to improve the ab initio results. 
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I. INTRODUCTION 

The study of the superheavy elements is now a popular 
area of research driven by the search for the island of 
stability in the region Z=\\A to Z = 126 where shell 
closures are predicted (see, e.g., [l[). Elements up to Z = 
118 have been synthesized (see, e.g., @, Q and evidence 
for the naturally occurring element Z = 122 has been 
reported [4j. 

Experimental efforts are underway to measure the 
spectra and chemical properties of the superheavy el- 
ements Jj|. There are also many theoretical works in 
atomic physics and quantum chemistry devoted to these 
studies (see, e.g., Refs. in@, Q). 

In our previous work [8j we calculated the electronic 
spectra of the superheavy element Z = 119 and the 
singly-ionized superheavy element Z = 120. The nu- 
cleus with Z = 120 protons and N ~ 172 neutrons is 
predicted to be doubly magic in relativistic mean-field 
nuclear calculations Q . Moreover, there is evidence from 
fusion reactions for the enhanced stability of the element 
Z = 120[ljl. 

To the best of our knowledge the spectrum of the neu- 
tral element Z = 120 has never been calculated before. 
On the other hand, it has a relatively simple electronic 
structure similar to the structure of its lighter analogues 
barium and radium. These atoms have two valence elec- 
trons above closed shells. The calculations for the atoms 
are possible with very high accuracy [ill [l2j by com- 
bining the correlation potential method with many-body 
perturbation theory and the configuration interaction 
technique. 

In the present work we perform accurate relativistic 
calculations for the energy levels of the neutral super- 
heavy element 120 applying a similar approach as in ear- 
lier works for barium and radium [111 LL31 . 



II. METHOD 

Calculations are performed with a method that com- 
bines the configuration interaction (CI) technique with 
many-body perturbation theory (MBPT). It has been de- 
scribed in detail in our previous works [H mill Qui. 



Here we repeat the main points, focusing on the details 
specific to the current calculations. 

Calculations are carried out in the V N ~ 2 approxima- 
tion This means that the initial Hartree-Fock pro- 
cedure is performed for the doubly-ionized ion, with the 
two valence electrons removed. This approach has many 
advantages. It simplifies the inclusion of the core- valence 
correlations by avoiding the so-called subtraction dia- 
grams [lH [HI . This in turn allows one to go beyond sec- 
ond order in many-body perturbation theory in the treat- 
ment of core- valence correlations. Inclusion of the higher- 
order core- valence correlations significantly improves the 
accuracy of the results 0, [IB] . 

The effective CI Hamiltonian for an atom with two 
valence electrons is the sum of the two single-electron 
Hamiltonians and an operator representing the interac- 
tion between the valence electrons, 



H cS = h x {r x ) + h x {r 2 ) + h 2 { ri ,r 2 ) 



(1) 



The single-electron Hamiltonian for a valence electron 
has the form 

hi = h + Si , (2) 
where ho is the relativistic Hartree-Fock Hamiltonian, 
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col ■ p + (/? — l)mc 2 h V 



N-2 



(3) 



and Si is the correlation potential operator which rep- 
resents the correlation interaction of a valence electron 
with the core. 

The interaction between valence electrons is given by 
the sum of the Coulomb interaction and the correlation 
correction operator E2, 



h 2 = 



ri - r 2 



£ 2 (ri,r 2 ) 



(4) 



The operator £2 represents screening of the Coulomb in- 
teraction between valence electrons by core electrons. 

The two-electron wave function for the valence elec- 
trons \t can be expressed as an expansion over single- 
determinant wave functions, 



(5) 
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TABLE I: Screening factors fk for inclusion of higher-order 
correlations into the exchange part of Ei and into £2 as func- 
tions of the multipolarity k of the Coulomb interaction. 



k 





1 


2 


3 


4 


5 


^cxch 

£ 2 


0.72 
0.90 


0.62 
0.72 


0.83 
0.98 


0.89 
1.00 


0.94 
1.02 


1.00 
1.02 



The functions <&i are constructed from the single-electron 
valence basis states calculated in the V N ~ 2 potential, 

$i(n,r 2 ) = -y=(ip a (ri)ipb(r2) - ^(nMiCo)) . (6) 

The coefficients Ci and two-electron energies are found by 
solving the matrix eigenvalue problem 

(H cS - E)X = , (7) 

where Hff = (QAH^i) and X = {c 1; c 2 , . . . , c„}. 

The most complicated part of the calculations is calcu- 
lation of the correlation correction operators Si and £2. 
We use MBPT and the Feynman diagram technique to 
do the calculations. The MBPT expansion for £ starts 

from the second order in the Coulomb interaction. In- 

(2) (2) 

elusion of the second-order operators £^ ' and £ 2 m t° 
the effective Hamiltonian (|T|) accounts for most of the 
core- valence correlations. However, further improvement 
is achieved if higher-order correlations are included into 
£1 and £2. 

We include higher-orders into £1 in the same way as 
for a single valence electron atom [l6|. Two dominating 
classes of higher-order diagrams are included by applying 
the Feynman diagram technique to the part of £1 that 
corresponds to the direct Coulomb interaction. These 
two classes correspond to (a) screening of the Coulomb 
interaction between valence and core electrons by other 
core electrons and (b) the interaction between an electron 
excited from the core and the hole in the core created by 
this excitation [l6[ 

The effect of screening of Coulomb interaction by the 
core electrons in the exchange diagrams is approximated 
by introducing screening factors fk (see Table UJ into 
each line describing the Coulomb interaction between 
electrons in a Brueckner-Goldstone diagram. We assume 
that screening factors fk depend only on the multipolar- 
ity of the Coulomb interaction k. The screening factors 
were calculated in our early works [H, [I?) and then used 
in a number of later works. It turns out that screening 
factors have very close values for atoms with similar elec- 
tron structure. In particular, the same values can be used 
for all atoms of the first and second columns of the pe- 
riodic table. The screening factors for £° xch were found 
by calculating the direct part of £1 with and without 
screening. 

A similar way of approximate inclusion of higher-order 
correlations via screening factors was used for £2. The 
values of the factors, however, are different (see Table HJ. 



These factors were found by comparing £1 in second- 
order and in all-orders with both screening and hole- 
particle interaction included. 

One needs a complete set of single-electron states to 
calculate £ and to construct the two-electron basis states 
([5]) for the CI calculations. We use the same basis in both 
cases. It is constructed using the £?-spline technique [HI, 
Ojj . We use 40 _B-splines of order 9 in a cavity of radius 
i? max = 40as, where as is Bohr's radius. The upper and 
lower radial components R^ ,l (r) of the Dirac spinors for 
single-electron basis orbitals ip a in each partial wave are 
constructed as linear combinations of 40 i?-splines, 

40 

R u J(r) = Y,l> U jBi{r) . (8) 

i=l 

The coefficients b^f are found from the condition that tp a 
is an eigenstate of the Hartree-Fock Hamiltonian ho ([3]). 

III. RESULTS 

Our results for barium, radium, and element 120 are 
presented in Table [Hi Ab initio results are listed in the 
column under the heading "£" . Here, the all-orders cor- 
relation potential £1 is used. This is the same £1 used in 
our previous work for Ba + , Ra + , and the singly ionized 
element 120 [8j. For both barium and radium we see that 
all ab initio calculated levels lie deeper than experiment. 
That is, the removal energies for the two s electrons from 
the ground states 6s 2 for barium and 7s 2 for radium are 
larger than the experimental results, both by 0.2%. The 
ab initio excitation energies are smaller than experiment. 
For barium the disagreement with measured values is on 
the order of 1%; deviations for the 6s5d states are larger 
than for the 6s6p states. The results for radium are bet- 
ter. 

As a way to empirically correct the results for radium 
and element 120, we use fitting factors found by fitting to 
measured levels for barium and radium, respectively. The 
empirical fitting factors /e a and /R a are different for the 
different partial waves s, p, and d and placed before the 
associated £1. This is done in the same way as in our 
work for single- valence electron systems Q. However, 
here the fitting is performed for the two valence elec- 
tron atom; it means that, as well as correcting for unac- 
counted core- valence correlations, information about the 
valence- valence correlations, incompleteness of the basis, 
and other higher-order corrections are also contained in 
this fitting. 

Results for spectra calculated with the empirical fit- 
ting factors are presented in Table [Til under the head- 
ings "/B a £" and "/R a £"- We see that by using the 
three fitting factors, energy levels for all states consid- 
ered for both barium and radium are in excellent agree- 
ment with experiment. This includes the singlet states 
Pi and 1 I?2- In our previous works on barium and ra- 
dium [ill [l2| , where similar methods were employed, the 
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TABLE II: Energy for removal of two s electrons (for So, in atomic units) and excitation energies (in units cm - ) for barium, 
radium, and element 120. Results of calculations with all-orders correlation potential appear in the column "E"; those with 
empirical fitting factors in the columns "/Ba£" and "/rtaE". ^-factors for element 120 are presented in the last column. 
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Ref. . 



singlet states were much worse. This improvement comes 
about from the different screening factors used in S2 in 
the current work. 

Now we turn to the results for radium and element 120 
found with these fitting factors. It is seen from Table 
HT1 that with empirical fitting factors from barium, all 
levels for radium move in the direction of experiment. 
However, the corrections are over-shot. The triplet 7s7p 
levels are seen to improve, from about -1% deviation to ~ 
0.1%. On the other hand, the 7s6d levels get worse, with 
deviations about 1%. We note that the largest deviations 
from experiment for the ab initio results for barium are 
for the 6s5d levels, with deviations —(3-4)%. The fitting 
factor for the <i-wave has clearly over-compensated these 
corrections. 

For element 120, the removal energy for the two s elec- 
trons from the ground state 8s 2 decreases by 0.2% with 
fitting. This is the same difference we saw between ab 
initio and measured values for barium and radium. For 



the excitation energies, the differences range from 0.1% 
for the higher levels to over 1% for the two lowest levels, 
with all fitted values larger than the ab initio ones. We 
see that the differences for element 120 are smaller than 
for radium. This is because the ab initio values for ra- 
dium are very good, better than those for barium. We 
take the fitted results as our final values for the low-lying 
spectra of element 120. 

The differences between the ab initio and fitted values 
gives an indication of the accuracy of the results. While 
for the higher levels these differences are about 0.1%, 
from a consideration of the large differences we saw for 
radium, we will assign a conservative uncertainty of 1% 
to our results for excitation energies. 

We considered Breit and radiative corrections in our 
work on the spectra of superheavy elements 119 and 120 + 
Q . We found that by fitting the results of pure correlated 
calculations to the experimental values for Cs and Ba + , 
the extrapolated values for Fr and Ra + agreed with ex- 
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periment at a level better than the estimated size of Breit 
and radiative corrections. This method was then used for 
the heavier homologues element 119 and 120 + as a way 
of taking into account these corrections semi-empirically. 

In the present work we also fit pure correlated values 
of the lighter homologues to experiment for use in calcu- 
lations for the heavier elements. In principal, this may 
be a way of taking into account Breit and radiative cor- 
rections. However, in this case our calculations are not 
at the same level of precision as in our work Ref. Q 
dealing with one-valence electon atoms and ions. In the 
current work, the uncertainty in the final results stem- 
ming from correlations is larger than the expected size 
of Breit and radiative corrections. Therefore, we do not 
need to account for these corrections. 

The ordering of the levels for barium, radium, and el- 
ement 120 are different. Relativistic effects increase ap- 
proximately as Z 2 a 2 and are larger for the heavier ele- 
ments. They are responsible for pulling in s and p levels 
which effectively screen d levels, pushing them out. This 
is the reason that 7 's7p 3 Pq is the first excited state rather 
than 7s6d 3 Di for radium. For element 120, stronger rel- 
ativistic effects are responsible for pulling the 8s8p 3 Pi 
level below 8s7d 3 D 1 . 

Note that we use the LS notation for the states. While 
the atoms under consideration are highly relativistic and 
the jj notation is more appropriate, the LS notation is 
used extensively in the literature (including in the tables 
of Moore [2(|) and for this reason we adopt it here for 
easier comparison. In this work we calculated g-factors 
for element 120 to help in identification of the states; 
these are presented in the final column in Table HT1 

We found in our work Ref. Q that there is a sizable 
volume isotope shift for the elements 119 and 120 + . In 
this work, as in Ref. calculations have been per- 

formed using a two-parameter Fermi distribution for the 
nuclear density with half-density radius c = 8.0 fm and 
10-90% width t = 2.0 fm corresponding to a root-mean- 



square charge radius r rms « 6.42 fm. We define the vol- 
ume isotope shift as 



where here the energies E are taken relative to the con- 
tinuum (removal energies for two electrons). Values of k 
for element 120 are listed in Table ITTT1 A table of values 
for r rms calculated in the nuclear Hartree-Fock-BCS ap- 
proximation yields a range from around r rms = 6.45 fm 
to r rms = 6.95 fm for the light to very heavy isotopes 
[2ll ]. For r rms = 6.90 fm, for example, the value for the 
removal energy for the two s electrons in the ground state 
is 8.5 x 10~ 4 a.u. less than for r rms = 6.42 fm. This differ- 
ence amounts to 0.1%. For the level 8s8p 3 Pq, the shift is 
5.0 x 10" 4 a.u. = 110cm -1 which is ~ 0.5% of the value 
for the excitation energy. 

TABLE III: Volume isotope shifts k for states of element 120. 

State 8s 21 S 8s8p s P 8s8p 3 Pi 8s8p 3 P 2 8sSp 1 P 1 
k -0.018 -0.012 -0.012 -0.012 -0.013 



IV. CONCLUSION 

The energy levels of low states of the superheavy ele- 
ment Z = 120 have been calculated with an uncertainty 
of about 1%. These results may be useful for experimen- 
tal studies of this element and for predicting its chemical 
properties. 
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